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Abstract 


Given a polynomial 


f{x) = aox"+aix" -h Un 


with positive coefficients at, and a positive integer M < n, we define a(n infinite) generalized Hurwitz 
matrix HM{f) ■ ■ We prove that the polynomial f{z) does not vanish in the sector 


^zeC: |arg( 2 )| < 

whenever the matrix Hm is totally nonnegative. This result generalizes the classical Hurwitz’ Theo¬ 
rem on stable polynomials (M = 2), the Aissen-Edrei-Schoenberg-Whitney theorem on polynomials 
with negative real roots (M = 1), and the Cowling-Thron theorem (M — n). In this connection, we 
also develop a generalization of the classical Euclidean algorithm, of independent interest per se. 


Introduction 

The problem of determining the number of zeros of a polynomial in a given region of the complex plane 
is very classical and goes back to Descartes, Gauss, Cauchy [3], Routh [221123) . Hermite [5], Hurwitz 
and many others. The entire second volume of the delightful Problems and Theorems in Analysis by 
Polya and Szegd m is devoted to this and related problems. See also comprehensive monographs of 
Marden [T^, Obreshkoff [19], and Fisk [6]. 

One particularly famous late-19th-century result, which also has numerous applications, is the Routh- 
Hurwitz criterion of stability. Recall that a polynomial is called stable if all its zeros lie in the open left 
half-plane of the complex plane. The Routh-Hurwitz criterion asserts the following: 

Theorem 1 (Routh-Hurwitz [T31[m[53]). A real polynomial f{x) = -|-oia;"“^ H-h a„ (oq > 0) 

is stable if and only if all leading principal minors of its Hurwitz matrix i? 2 (/) up to order n are positive. 

Decades after Routh-Hurwitz, Asner |5| and Kemperman m independently realized that the Routh- 
Hurwitz criterion can be restated in terms of the total nonnegativity of the Hurwitz matrix. Moreover, 
the Hurwitz matrix of a stable polynomial admits a simple factorization into totally nonnegative factors 
m- These developments are described in dSl Section 4.11]; see also a separate section j^Hl Section 4.8] 
on generalized Hurwitz matrices. The converse direction of the total nonnegativity criterion was fully 
established only a few years ago in m- 

Theorem 2 ([IS])- A polynomial f{x) = aox" -I- aix" ^ + • • • + Un {clq, oi,..., a„ S K; cq > 0) has no 
zeros in the open right half-plane Rez > 0 if and only if its Hurwitz matrix i? 2 (/) is totally nonnegative. 
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Note that the strict stability has to be replaced here by its natural weaker counterpart. As we shall 
discuss in Section [SI generalized Hurwitz matrices turn out to enjoy analogous properties! (We should 
add that Pinkus EOl Section 4.11] mistakenly asserts that the total nonnegativity of the infinite Hurwitz 
matrix is not sufficient for all zeros to lie in the closed left half-plane. However, it is the finite, i.e., nxn, 
Hurwitz matrix, whose total nonnegativity is not sufficient.) 

Another famous mid-20th-century result of Aissen, Edrei, Schoenberg and Whitney [I] concerns a 
seemingly different class of polynomials, i.e., those with real negative roots. This result admits a strikingly 
similar formulation to the total nonnegativity Theorem |2] 

Theorem 3 ([I]). A polynomial f{x) = oox" -I- H- -f a„ (oq, oi,..., a„ S K; oq > 0) has only 

real negative zeros if and only if its Toeplitz matrix Hi(/) is totally nonnegative. 

The reader may wonder why we refer to the Hurwitz and Toeplitz matrices associated to a polynomial 
/ as H 2 {f) and Hi{f). As we shall see, both are special cases of generalized Hurwitz matrices, which will 
be denoted by HM{f), for M = 1 and M = 2. Note that we take all matrices Hm to be infinite. Fittingly, 
our main result generalizes both Routh-Hurwitz and Aissen-Edrei-Schoenberg-Whitney criteria: 

Theorem 4. A polynomial f{x) = qqX^ + aix^~^ -I- • • • -f a„ (oq, oi,..., On £ R, oq > 0) of degree n 
has no zeros in the sector 

{zGC: |arg(z)| < 

whenever its generalized Hurwitz matrix Hm^J) is totally nonnegative. 

While not enjoying the ideal ’if and only if’ format, this result is nevertheless beautifully similar to 
the two classical results we just revisited. We shall discuss obstacles to the converse statement for M > 2 
in Section [51 Section |T] is devoted to a generalized Euclidean algorithm, which is crucial to our proof. 
That algorithm itself should be of independent interest in algebra and analysis. In Section |21 we shall 
examine various continued fractions that can be constructed from the generalized Euclidean algorithm. 
In Section [3] we shall properly introduce generalized Hurwitz matrices and point out their connections 
with the generalized Euclidean algorithm. We shall then examine the interplay between generalized 
Hurwitz matrices and their regular Hurwitz submatrices in Section j?) The proof of our main theorem 
will be given in Section [51 along with an interesting factorization of generalized Hurwitz matrices in 
Section ini Various related questions will be discussed in the remaining two sections. 

Our approach is based on a systematic study of the following objects: 

• generalized Hurwitz matrices and their submatrices; 

• continued fractions; 

• roots and coefficients of polynomials, especially Routh’s (or Euclidean) algorithm. 

We shall begin by defining and examining the generalized Euclidean algorithm, which will play a crucial 
role connecting all objects mentioned above. 

1 Generalized Euclidean algorithm 

We will now develop a generalization of Euclidean algorithm for M > 2 polynomials (replacing M = 2 
for the regular Euclidean algorithm). What folows is in fact a special version designed for the purposes 
of splitting a given polynomial into M parts according to the residues of the coefficients mod M. A more 
general version of the generalized Euclidean algorithm will be described in another paper. 

Let 

f{x) = qqx'^ -I- aix’^~^ + On 

be a polynomial of degree n with real coefficients aj, j = 0,... ,n. As usual, we define deg f{x) : = — oo 
if f{x) = 0. Let M be a positive integer, 2 < M < n. Then / can be split into a sum of polynomials 

fix) = foix) + fi{x)-\ - + fM-iix), (I) 
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where 


( 2 ) 


fji^) = 

l=j mod(Af) 

0<J<n 

Definition 5. A polynomial p{x) is called arithmetic with difference M (and residue k) if 

p{x) = ^ 

l=k m.od(M) 

Per this definition, (ED is a decomposition of / into a sum of arithmetic polynomials with difference M 
and the full set of residues. The degrees of monomials in every fj{x) form an arithmetic progression with 
difference M; every nonzero monomial term of f{x) enters only one fj{x). Notice that a zero polynomial 
is arithmetic for any difference M. 

A generalized Euclidean algorithm associated with the decomposition ([T|) is defined as follows. For 
any i = 0,1,..., M — 2, there is a unique representation of the form 

fi{x) = di{x)fi+i[x) + fi+Aiix), 

where di(x) and fi+M{x) are polynomials, subject to the following rules: 

(a) If deg(/i) > deg(/i+i) > —oo, then fi+M is the remainder in the division of the polynomial fi{x) 
by fi+i{x), and di{x) is the quotient. Hence deg(/i+M) < deg(/i+i). 

(b) If -oo < deg(/i) < deg(/i+i) then di{x) = 0 and fi+M{x) = fi{x). 

(c) If -00 = deg(/i+i) then di{x) = 0 and fi+M{x) = Mx). 

The algorithm stops when fn(x) is constructed. 

We call this process the generalized Euclidean algorithm with step M applied to the polynomials 
{/o, /i, • ■ •, /m-i}- 


Example 6. Let M = 3 and f{x) = x'^ + x^ -\- x^ + x'^ + + x'^ + x + 1. Then 

/o(x) = + x'* + x; 

/i(x) = X® + x^ + 1; 

/2(x) = X® + x^. 


Next, 


/o(x) = X - /i(x) +0 
/l(x) = X • / 2 (x) + 1 
/ 2 (x) =0-0 + / 2 (x) 
/ 3 (x) =0-1 + 0 
fi{x) = 0 - / 5 (x) + 1 


h{.x) = 0 
/4(x) = 1 
/5(x) = X® +X^ 
h[x) = 0 
h{x) = 1. 


It is convenient to arrange the resulting polynomials into the following table: 


Groups: k 

0 

1 

2 

0 + 3fc 

x‘ + x"^ + X 

0 

0 

1 + 3fc 

X® + x^ + 1 

1 

1 

2 + 3fc 

X® + x^ 

X® + x^ 



Example 7. Let M = 3 and /(x) = (x + 1)’^ = x’^ + 7x® + 21x^ + 35x'^ + 35x^ + 21x^ + 7x + 1. Then 

/o(x) = x^ + 35x'^ + 7x; 

/i(x) = 7x® + 35x^ + 1; 

/ 2 (x) = 21x® + 21x^. 
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Next, 



= y • /i(a;) +30a;^ + yx 

fsix 

fi{x) 

= 1 • h{x) +28a;® + 1 

fi{x 

h{x) 


^ fdix 

h{x) 

= ^ ■ /*(») + S" 

h{x 

f^ix) 

5 • 28a: , , 

= •/5(a;) + l 

h{x 


The table for this polynomial looks as follows: 


30x‘^ + 
28x^ + 1 


81 x^ 

5 

81 

14^ 

1 . 


Groups: k 

0 

1 

2 

0 + 3fc 

1 + 3fc 

2 + 3fc 

x^ + 35a; + 7x 
7x® + 35a;® + 1 

21a;® + 21a;2 

30a;4 + fx 
28a;® + 1 

Sla;^ 

_5_ 

rp 

14 

1 


Example 8. Let M = 3 and f(x) = x'^ + x^ + x^. Then 

fo(x) = 
fi(x) = x®; 
/2(x) = x®. 

Next, 


fo{x) = 

= X ■ 

' fi{x) 

+ 0 

h{x) 

= 0 

fl{x) = 

= X ■ 

' f2{x) 

+ 0 

f4{x) 

= 0 

f2(x) = 

= 0- 

faix) 

+ a® = 

^ f5{x) 

= x® 

fsix) = 

= 0- 

h{x) 

+ 0 

h{x) 

= 0 

f4(x) = 

= 0- 

f5{x) 

+ 0 

h{x) 

= 0. 


The table for this example is the following: 


Groups: k 

0 

1 

2 

0 + 3k 

x'’ 

0 

0 

l + 3k 

X® 

0 

0 

2 +3k 

X® 

X® 



Example 9. Let M = 4 and f{x) = x^ + x'^ + x^ + x^ + x^ + x^ + x. Then the table looks like 


Groups: k 

0 

1 

2 

0 + 4:k 
l+4:k 

2 + Ak 

3 + 4A; 

X^ + X® + X 

0 

X^ + X® 

X® + x^ 

x'® + X® + X 

0 

0 

X® + x^ 

(1 c 

X + X + X 

0 


Theorem 10. Let f{x) be a polynomial of degree n > 2 with real coefficients and let 2 < M < n. Let 
{fo, fi, ■ ■ ■, fn} be the polynomials obtained by the generalized Euclidean algorithm with step M. Then 

A. Every nonzero fi{x) is arithmetic with difference M. 

B. If both polynomials fi{x) and fi+M(x) are nonzero, then their residues are equal. 
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C. For every i either di{x) = 0 or di(x) is an arithmetic polynomial whose residue equals the difference 
between the residues of fi and fi+i- In particular its degree is greater or equal to that difference. 

D. If fi{x) and fi+i{x) are nonzero polynomials, then 


|deg(/i)-deg(/i+i)| > 1. (3) 

Proof. Consider the table of f{x) associated with the generalized Euclidian algorithm with step M: 


Groups: k 

k = 0 

k = I 

k = 2 


k =last 

0 + Mk 

1 -h Mk 

fo(x) 

fi(x) 

fnix) 

fi-\-M{x) 

f2M{x) 

fl+2M{x) 


fn-jix) 

fn-j+1 (x) 

j -\- Mk 

fjix) 

fj-i-Atix) 

fj+2Mix) 


fn(x) 

M{k + l)-l 

fM-i{x) 

fM-l + M{x) 

fM-l+2Mix) 




By the definition of the generalized Euclidean algorithm, all nonzero fi{x) with 0 < * < M — 1 in 
group 0 are arithmetic with difference M and distinct residues. It follows that no pair of them can have 
equal degrees, implying (|31) ioi i = 0,... ,i = M — 2. 

Group 1 of the polynomials fi+M, i = 0,. ■ ■, M — 1 (see column I in the table above) is determined 
by the identities 

fi{x) = d^{x)f^+i{x) + fi+M{x), i = 0,1,... ,M - 1. 

IfO<i<M — 2 and we have the case (a) of the algorithm, then fi+M is the remainder in the 
division of one arithmetic polynomial, fi{x), by another arithmetic polynomial, fi+i{x) (with a different 
residue). It follows that di{x) is an arithmetic polynomial with residue k > 0, which shifts the arithmetic 
progression of the exponents in fi+i(x) to the arithmetic progression of the exponents of fi{x). Hence 
fi+M = fiix) — di{x)fi+i{x) is either the zero polynomial or an arithmetic polynomial with the same 
residue as fi{x) and deg(/i+M) < deg(/i+i) < deg(/i). 

In cases (b) and (c) we have fi+M = fi- Moreover, if fi{x) = 0, then by requirement (b) of the 
algorithm all polynomials fi+kM, k = 1 ,are zero too. If fi+i{x) = 0, then by requirement (c) of the 
algorithm fi = fi+M- Since fi+M+i = 0 by our argument above, we get fi = fi+M = fi+ 2 M- Continuing 
by induction, we see that the whole row starting with fi is filled with fi. 

As to the ’boundary’ pair /m-i and /m, in case both of them are nonzero, the residue of /m equals 
the residue of /o, which is n, whereas the residue of /m-i is n — M + I. Since n — {n — M + 1) = M — 1, 
we see that polynomials /m-i and /m have different residues, implying Q in this case too. 

The process can be continued by induction on k. In this way, we obtain two possibilities for each 
row: In the first case we have polynomials of strictly decreasing degrees and the row terminates either 
in a string of zero polynomials or in several copies of the same polynomial. In the second case the row 
is simply filled with zeros. In the third case it is filled with copies of the same polynomial. 

Corollary 11. Let f{x) be a polynomial of degree n > 2 with real coefficients and let 2 < M < n. 
Let {/o, /i,..., fn} be the polynomials obtained by the generalized Euclidean algorithm with step M. 
If none of these polynomials is zero, then di{x) = CiX for i = 0,1,..., n — 1 and deg(/fe) = n — k. 

Remark. We will refer to the situation of Corollary ITT] where none of the polynomials {/o, /i, • ■ • /«} is 
zero as the non-degenerate case of the generalized Euclidean algorithm. 


2 Continued fraction expansions 

Of course, the regular Euclidean algorithm can be also applied to any pair of polynomials (fi, fj) gener¬ 
ated by the generalized Euclidean algorithm. Let us look into this, assuming the non-degenerate case. 
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fl 

Suppose 0 < i < j < M. Denote the fraction ^ by Rij. We want to represent the function Rij{z) 


as a continued fraction 




= qf{z) + 




( 4 ) 


+ 




Applying the ordinary Euclidean algorithm to the pair (/j, we construct a sequence of polynomials 
fo^ fW ■ • ■) fk leading coefficients h^, , ..., , respectively, as follows: 


foix) ■■= fiix)=a,x'^ * + ai+Ma;” * ^ + a^+2Mx'^ * 


fmx) := fj{x)=ajX^^+aj+MX 


n—j — M 


aj+2MX 


,n—j — 2M 


The subsequent polynomials are defined by extracting the leading term from the ratios 


cv 

Ji-i 




= 1 , 2 ,... 


where 


gfix) := 


' /jb 

jj: £ jg even; 

hV 


- x-^ 

hV 


if £ is odd. 


degfil^ix) =n-j-M 
Hence the continued fraction (jT|) can be written explicitly as 


i?y(z) = 




!h_ , 


"■k-l 

"-k 


zH 


where 


M — {j — i) if k is even; 


= 


J - * 


if k is odd. 


( 5 ) 


We now proceed to make a simple but crucial observation about the continued fractions Rij. It turns 
out that these fractions, viewed as functions on C, map cones with sufficiently small apertures emanating 
from the origin to similar cones. Here are the details. 

Given two angles a < /3 between — tt and tt, consider the cone 

Ka,i 3 '■= {z G C : a < arg(z) < /3}. 

Lemma 12. Let i? be a continued fraction 

R{z) = aiz^ + -^ 


a2Z 


M — m 


1 


1 


asz' 


QkZt^ 


with all coefficients aj positive, j = 1,..., fc, with m an integer between 0 and M, and with the exponents 
alternating between m and M — m (so that /r = m if fc is odd and = M — m ii k is even). Then, for 
any angle a G [0,7r/M], the function R maps the cone Ko a into the cone K_(^M-m)a,mcf 
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Proof. The map i? is a composition of special monomial maps, multiplication by constants, additions, 
and inversion. Let us examine how these maps act on our cone Kq ^- 

A function ()^ ■. z ^ maps a cone to Multiplication by a positive constant leaves 

any such cone invariant. Inversion ()“^ maps a cone Ka^p (of course excluding the origin) to the cone 
K-p^-a- Equipped with these basic observations, we can now understand the action of i? on a cone of 
type ATo.a for a < tt/M. 

Suppose for simplicity that k is odd. Then the last monomial. O’", maps ATo.a to Ko^rna- Multipli¬ 
cation by Ofc leaves the latter cone invariant, and inversion maps it to the cone AT-ma.o- The previous 
monomial (followed by muplication by a positive constant Uk-i) maps Kq ^ to so the result 

lies in KQ^(^M—m)a T ^ ^—maPM—m)a- 

The inversion that follows maps K_jy^a,{M-m)a to The function afc_ 2 ()’" maps ATg^a 

to ATo.ma, and the next cone addition yields C K-{M-m)a,ma- From this point 

on, we shall alternate between the cone and its reflection since 


KopM-m)a + (A1_(m — m)a,ma ) 


-1 


K, 


0,7710; 


+ iK_ 


ma,(M—m)a 


^ ^ —ma,{M—m)a 

^—{M—m)a,Tna- 


Our last operation is of the second type, so we shall end up inside the cone K_(^M-m)a,ma- Note that 
this happens regardless of the parity oi k. □ 

As promised above, we can now apply this lemma to our functions Rij : 


fi 

Corollary 13. Let Rij = -A (0 < * < j < M— l)bea rational function defined 

Jj 

fraction ® with all coefficients r = 0, ..., k positive. Then Rij maps the cone 
K_(^M-j+i)a,{j-i)a whenever 0 < a < tt/M. 


by continued 
into the cone 


3 Generalized Hurwitz matrices 


Every polynomial 


f{x) = oga;" -I- aix" ^ + ■■■+ an 


and an integer 2 < M < n determine a generalized Hurwitz matrix (defined by Goldman and Sun in [8]) 


Hm = {aMj-i)^j=l- 

In the formula for the entries of Hm we assume that = 0 for integers fc < 0 and k > n. Thus a 


generalized Hurwitz matrix is an infinite matrix of the following 

form: 


/ om-i 

a-2M-i 

asM-i 

...X 


aM-2 

a2M-2 

a^M-2 



ao 

O-M 

a2M 



0 

om-i 

a2M-i 


Hm = 

0 

aM-2 

a2M-2 



0 

Og 

Om 



0 

0 

a-M-i 



V : 



■■■/ 


The matrix Hm is constructed from the coefficients of the arithmetic polynomials /m-i, /m- 2 , ■ • /o 
of the polynomial /. The first row is filled with the coefficients of /m-i- The second row is filled with 
the coefficients of /m- 2 , etc. The coefficients of /g form the Mth row of Hm- Then this first block of 
M rows is shifted one step to the right and placed underneath. And so on. Note that this structure 
generalizes both Toeplitz and Hurwitz structures simultaneously. 
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Definition 14. A matrix A of size n x m, where n and m may take infinite values, is called totally 
nonnegative if all its minors are nonnegative: 


A 


*1 



> 0 . 


Taking p = 1, we see that all entries of a totally nonnegative matrix are nonnegative. 

It was shown in [5] that the generalized Hurwitz matrix Hm is totally nonnegative if and only if its 
n special minors are nonnegative: 


HM{k,r): =Hm 


fc + 1 
2 



( 6 ) 


for k = 1,M — 1 and r = 1,.. 


n+fc—1 I 

M-1 J * 


This enumeration of special minors was introduced by 


Pinkus in [20l Theorem 4.6, p.lll]. 

Another natural enumeration was introduced by Goldman and Sun in [5]. For this enumeration, we 
need to locate the bottom left matrix entry ap of the minor ([S]) (i = k r — 1 and j = r): 


p = Mj — i = {M — l)r — (fc — 1) = [M — l)(r — 1) + (M — k). 


(7) 


Every number p = 1, 2,..., n can be uniquely represented in the form ([7]). Therefore if we are given such 
a p, we can determine k and r. We denote by 


Ap:=HM{k,r) = Hm 


k + 1 
2 


fc + r-1) 


the corresponding minor of the generalized Hurwitz matrix. 

Theorem 15 ([8l Theorem 2.1]). Suppose that all special minors 

Ap > 0, p=l,...,n (8) 


are positive. Then 

(a) the matrix Hm is totally nonnegative; 

(b) all coefficients ap of the polynomial f{x) are positive; 

(c) a minor M of Hm is strictly positive if and only if its diagonal elements are strictly positive. 

Our next theorem shows how generalized Hurwitz matrices are related to the generalized Euclidean 
algorithm. 

Theorem 16. If all the leading coefficients {Hq, hi, ..., /i„} of the polynomials {/o, /i, ..., fn} are 
positive, then the conditions ([5]) are satisfied and the matrix Hm is totaly nonnegative. 

Corollary 17. If all the leading coefficients {ho, hi, ..., /i„} of the polynomials {/o, fi, ■ ■ ■, fn} are 
positive, then all coefficients of the polynomial / are positive. 

The proof of Theorem [TH] is an easy consequence of Theorem [TS] and the following lemma. 

Lemma 18. Let oq > 0 and let hi denote the leading coefficients of the polynomials fi in the generalized 
Euclidean algorithm with step M. Then 

= hM-lh2M-2 ■ ■ -hrM-r 
= ^M-2^2M-3 • • • ^rM-(r+l) 


HM{M—l,r)— hihM ■ ■ ■ h^M-{r+M-2)t 


where r = 1,..., 


M-1 


, hi 0 for i > n. 








Proof. The matrix Hm is made of the coefficients of polynomials /o, /i, •./m-i, arranged into 
shifted blocks: 


/ OM-I 


OpM-l 

...\ 

aM-2 

«2M-2 

a3M-2 


ai 

om+i 

a2M+l 


Op 

Om 

a2M 


V : 



■J 


By© there are iVf — 1 special minors of order r = 1. Their values are 


hn-i = om-i = 1) > 0, ..., hi = ai = Hm{M — 1,1) > 0, 

as is claimed by © for r = 1. 

There are M — 1 special minors of order 2 corresponding to the following matrices: 


1 a,M-i 

a2M-i\ 

f aM-2 

a2M-2 I 

\aM-2 

a2M-2 J ’ 

\aM-3 

a2M-3 J 


f Ol flM+i 
\ao om 


To evaluate the determinants of these matrices we apply Gauss elimination. Excluding ap in the last 
matrix using oi from the first row, we obtain a row equivalent matrix: 

(«! «M+l\ ^ «M+1 

\ao om J \ 0 Hm 


Indeed 

fo{x) = aox^ + H-; 

fi{x) = H-. 

Since fo{x) = ^xfi{x) + futix), we obtain that 


lilt \ n-M , ( \ n-2M , 

iMix) = 1 OM - —^M+1 j X + I a2M - —a2M+i jx H-. 

To summarize, this elimination results in annihilation of op, op —>■ 0, and in the replacement of all other 
coefficients in the row of op by the coefficients of the polynomial /m- In particular, om ^m- Similarly, 


f 02 

aM+2 \ 

( 02 

OM+2 \ 

1 Oi 

OM+1 ) 

AO 

hm+i J 


( OM-I 

02M-lI 

( OM-I 

02M-1 \ 

\a-M-2 

02M-2 / 

A 0 

h2M-2 J 


In other words, 0 and UM+k —^ hM+k- This proves © for r = 2. 

Let r > 2 and 1 < k < M — 1. The minor H{k, r) is the determinant of the square matrix 


1 aM-k 

02M-fc 

OsM-fe • 

• QrM—k 

Oi 

om+i 

02M+1 • 

0(r_l)M+l 

Op 

Om 

a2M 

0(r_i)M 

0 

om-i 

02M-1 • 

0(r_i)M-l 


V : : ■■ ) 

Eliminating op from the first column using oi, then ai using 02 , and finally aju-k-i using aM-k, we see 
that the above matrix is row equivalent to the matrix 


/aM-k 

02M-fc 

aSM-k 

arM-k \ 

0 

h2M-k-l 


h* 

'^rM-k-1 

0 

hM+1 

h* 

'^2M+1 

h* 

0 

Hm 

'^2M 

h* 

0 

OM-1 



V i 



/ 
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where the stars denote the coefficients of the polynomials f 2 M-k-i,- ■ ■, /m+i, /m- 

Observe that the second column of this matrix ends either with zeros or with some aj, M — 1 > j > 0. 
The rows of the matrix with the same indices make up the sequence of the coefficients of the polynomials 
fj, .. /m-1, • ■ •) f2M-k-i- It follows that we can run the same elimination process in column two as 
we did already in column one. As a result, we see that our matrix is row equivalent to 


initial matrix: 


/ QM-k 

a-2M-k 

asM-k 

QrM-k \ 

0 

h2M-k- 

h* 

1 H^M-k-l 

h* 

'^rM-k-1 

0 

0 

hzM-k-2 


0 

0 

^2M+1 


0 

0 

h2M 


0 

0 

h2M-l 


V : 




3 continued until 

we obtain a 

diagonal matrix w 

iM-k 

a2M-k 

asM-k 

* * ' OjrM—k 

0 

h2M-k-l 

h* 

h* 

'^rM-k-1 

0 

0 

hzM-k-2 


0 

0 

0 


0 

0 

0 


0 

0 

0 


0 

0 

0 

h-pM — k — {r—l 


It follows that 

which proves the lemma. 


HAfik, r) = hM-k ■ h2M-k-l ■ ■ ■ hrM-k-l, 


Corollary 19. Let oq > 0. Then the leading coefficients of the polynomials fi, i = 1, ..., n, satisfy 

HMik + l,r) 


where r = 


i 

M - 1 


k = r{M HM{k + l,0) := 1. 


Proof. Applying Lemma [T51 we obtain: 


h 


rM — r — 


h 


rM—(r+l) — 


HM{2,r-l) 


.rM-(r+M-2) (M - 1, r - 1) ’ 


Thus 


h 


rM—(r+fc) 


HM{k + l,r) 
HjAiik + l,r - 1)’ 


where 0 < k < M — 2. Assuming i = rM — (r + k), we obtain r 


i 

M - 1 


and k = r{M — 1) — i. 


□ 
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Example 20. Let us consider an example for the case n = 6, M = 3. Then 

f{x) = qqx^ + aix^ + a2x‘^ + a^x^ + 043 ;^ + a^x + ag. 
The generalized Hurwitz matrix looks like 

^02 05 0 0 .. A 

ai 04 0 0 ... 

= 0-0 0,3 06 0 ... 

0 02 05 0 ... 


There are six special minors: 

iL3(l,l)=a2, i^3(l,2) = 

iL 3 ( 2 ,l)=oi, H 3 ( 2 , 2 ) = 
By the generalized Euclidean algorithm, 

fo(x) = ooa;® + asx^ + oe 

fi(x) = aix^ + a^x^ 
f2ix) = a2x'^ + agx 


02 05 

Oi 04 


Oi 04 
Oo 03 


02 

05 

0 

^3(1, 3) = 01 

04 

0 

Oo 

03 

06 

Oi 

04 

0 

113(2, 3) = Oo 

03 

06 

0 

02 

O5 


Jz[x) = I 03- — j X +ae 

f4:[X) = 04- X 

V 02 / 


fdix) = 05- 


aia2tt6 

aias — aoa4 


We have: 


feix) = 06. 

ho = Oo, hi = ai = i?3(2,1), ft,2 = 02 = i?3(l, 1); 

I. 1 / ^ 1 oi 04 113(2,2) 

03 = —(0103 - 00O4 = — = TTTTTi’ 

oi oi Oo 03 113(2,1) 

j. 1 / ^ 1 02 05 ^^3(1,2) 

/14 = —(02O4 - O1O5) = — 

(I2 ^2 ^4 X 73 (l,l) 

, 11, ^ ^ 01 04 0 H3(2,3) 

05 = — — (050301 — 050004 — 010206) = —— Oo 03 06 = ,, ; 

01/i3 01/13 0 H 3 i 2 , 2 ) 


hg — ag — 


i?3(l,3) 


^3(1, 2)- 

4 Submatrices of generalized Hurwitz matrices 

Let HM(f) be a generalized Hurwitz matrix associated to a polynomial 

f(x) = uqx^ + aix^~^ H-+ o„. 

(a) 

We denote by its infinite submatrix determined by two polynomials fi and fj, 0 < i < j < M — 1: 


ttWI _ 
— 


/% 

OM+i 

G2M+j 

G3M+j --A 

Oi 

OM+i 

02M+i 

03M+i 

0 

aj 

GM+j 

02M+j 

0 

V; 


OM+i 

G2M+i 
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This matrix is the (ordinary) Hurwitz matrix of the polynomial 




A:>0 


Its arithmetic polynomials for M — 2 are given by 


prH^)=Y: 




%-2k r^iij) 




■,-2k-l 


k>0 


k>0 


Lemma 21. The even and odd parts of the polynomial satisfy 

(x^) = (x^) , pf^) (x^) = ^ 

Proof. By ([5]) and (flUl) for Pq*'^^ (x^) we have 


_ ^m-M-2{n—i) 




,2n —2(fcM+z) ^m-M—2{n—i) 


fiix )■ 


k>0 


Similarly, 




k>0 


( 10 ) 


as claimed. 


□ 


Lemma 22. Let f{x) = fo{x) + /i(x) + • • • + fM-i{x) and let all the coefficients Hq, ..., in the 
generalized Euclidean algorithm with step M be positive. The degree m of the polynomial P^‘^^\x) 

(a) 

whose Hurwitz matrix is is given by 



n — i 


if 

n — i 


r« - j] 

[ M \ 


[ M \ 


M 

b 

n — i 

+ 1 

if 

n — i 


n-j\ 

l M \ 

[ M J 


L M J 


(a) 

Proof. To find a formula for m we observe that is an infinite-dimensional submatrix of Hm 

which is obtained from Hm by keeping only rows corresponding to the coefficients of fi and fj. Thus, 
the principal minors of coincide with some minors of Hm- By Theorem llhl lcl. the leading principal 
minors of H^'^ are strictly positive if and only if their diagonal elements do not vanish. The diagonal 
elements of H^^ are: 


Gj, GjVf+2 5 ^2M-\-it ^2M-\-jt ^3Af+2 5 ^ZM-\-jt • • • • 


This sequence of positive numbers terminates as soon as either kM -|- i or kM j becomes greater than 
n for the first time. The entry akM+i is the last nonzero element if and only if 


kM < n < kM + j 


n-j 

M 


< k < 



n — i 


'n-j' 

M _ 


M 


since the difference between the bounds for k is strictly smaller than 1. The entry akM+j is the last 
nonzero element if and only if kM j < n < (k + 1)M -|- i or, equivalently, 


k < 


n-j 

M 


< 


n — i 

'W 


< fc -I- 1 


k = 


n — i 


n-j 

M 


M 


There are 2[(n — i)/M\ nonzero leading principal minors of H^P in the first case and 2[(n — i)/M\ -\-1 
nonzero leading principal minors in the second. By Routh-Hurwitz Theorem [TJ the number of nonzero 
leading principal minors equals the degree m of the polynomial [x ). □ 
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For our next (somewhat technical) lemma, we need to define a and /3 by 


M 


M 


, P = 


n-J 

M 


n-j 

M 




Lemma 23. Let £(m) = ^ 

3(b) 


■. Then 


PW; ^ ^e(m)-2ay^ ^ pfo) (j.) = j. . ^-e(m)-2p ^ 


Proof. By Lemma [21] 


By Lemma im 




m — 2- 




-2a 


if m is even 


M 1— 2a + l if m is odd, 
which proves the hrst formula. To prove the second one, we observe that 

'2 (|_^J - - 1 if w is even 


™ ~ 2 ——- 1 = 

M 


By Lemma [ 

n- ] 

I M \ M 

if m is odd. For even m we have 


([vj-v) if mis odd. 


= 2 


n-j 


n-j 


I M i M 


P = 


n — i 
M 


n- J 
M 


n-j 

M 


n-j 

M 


= -2/3. 


e(o,i). 


( 11 ) 


It follows that 


n-j 

M 


cannot be an integer and therefore 


'n-j' 


n-j 

M 


M 


implying the second formula. □ 

Corollary 24. Let f{x) = /o(x) + /i(x) + • • • + fM-i{x) and let all the coefficients ho, ..., h„ in the 
generalized Euclidean algorithm with step M be positive. Then 


(a) every polynomial 


p(b)(^) = P^^^\x) + P[^^\x) = (x^/^) 





is stable; 

(6) all coefficients in the (ordinary) Euclidean algorithm applied to the pair of polyno¬ 

mials {fi, fj) are positive. 

Proof. Item (a) follows from Lemma and the classical Routh-Hurwitz Theorem [T] 

To prove (b) we apply Lemma HI] with x replaced by y/x, x > 0. We have 


Rij (x) 


fj{x) p(i3)(^^M/2y 
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It is well known (see [THl Corollary 7.33]) that the C-fraction expansion of the rational function 




= diz 


d2Z+- 


d^z 


1 

d.z 


has all positive coefficients di, . ■ds whenever P{z) is a stable polynomial. It follows that 


Rij (x) — 


_ U-i)-M/2 _ Pq _ (j_i)_M/2 

(y). 


fjix) 

= + 


pWi(^M/ 2 ) 


= x^- 


{dicr^/2 


^22;M/2 


1 


= diX<-^-^^ + 




1 


I 


d^x^/"^ 


^32;M/2 

1 


d2X^-(^-^'^ + 


^M/2-ij-i) 

d^x^/"^ _l- 


d2X^~^^~'^'> 




Comparing this continued fraction with the continued fraction in ([5]) and invoking the uniqueness theorem 
for C-fractions (see m), we obtain that 


0 < di = 


h\^ 


-1 


hi 


hi > 0 , 


since = Ui > 0 and hi = aj >0. 


□ 


5 Main Results 


This section is devoted to the proofs of Theorem \T7\ for M > 2 and of our main Theorem Uj Note that 
these theorems differ both in assumptions and in conclusions. In particular, the statement of Theorem|T| 
includes the important case M = 1, which has to be excluded in the formulation of Theorem [571 We 
hope this does not unduly confuse the reader. We shall begin with Theoremand shall build the proof 
lemma by lemma, culminating in its proof. We shall then prove Theorem |4| using a different method. 

Given M >2 nonzero numbers w, zi, ..., zm-i in the cone = |u : 0 < arg(M) < —| and 

positive numbers a„ > 0, n > M, we define recursively 

Zn--=anWA - ^ -. ( 12 ) 

2'ri —1 • . . . • 

For n > M we denote 

S{n) : ={n — {M — 2),..., n} 

and define the partition S{n) = S-{n) U <S'+(u) of this set by 


S+{n) = {ie S{n) : arg(zi) > 0}, 
S-{n) = {i € S{n) : arg(zi) < 0}. 


Lemma 25. Let {zi, ..., z„} be a sequence defined by (fT2ll . Suppose that 0 < aTg{w) < and 


- TT + — < arg(z,) < —, i = M,...,n. 


(13) 


Then 


max 

S=S+{n),S-{n 


arg(zi) 

ies 


7 r(M — 1) 

< —^^ < TT. 

M 
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Proof. We use induction on n. Let n = M. Then S{M) = Since arg(zi) > 0 for 

i — — 1, we have either S+{M) = S{M) or S+{M) = S{M) \ {M}, hence S-{M) = 0 or 

S-{M) = {M}. Suppose hrst that S-{M) = 0. Since we assume that —tt < a,rg{zM) < tt/M, we obtain 

0 < arg(zM) < 


which implies that 


M 


X—' / \ X—' / \ TriAd — 1) 

2^ arg(zi) = 2^ arg(zi) < - tt- < tt. 


M 


ieS+(M) i=2 

Let us consider the case S-{M) = {M}. Then S+{M) = {2,..., M — 1}, implying that 


M-1 


o< i: «gfe)=E “«(*-•) 


ies+(M) 


M 


M 


For S-{M) = {M} we set 


ZM-l ■ ■■■■ Zi 


Then 


M-l 


(M — IItt g „ 

-arg(zi) = arg(M) < 0. 


Since S-{M) = {M}, we obtain 


Both w and u are in the half-plane 


arg(aMW + u) = arg(zM) < 0. 


(M — IItt tt 

-z :-TT-< arg(z) < — 

M - - M 


Therefore the sum qmw -I- m is in the same half-plane too. It follows that 

(M — l)7r 


M 


< SLTg{aMW + u) = arg(zM) < 0. 


This proves the Lemma for n = M. 

Suppose that the Lemma is proved for n. Let us prove it for n -I- 1. We have 


S{n) = {n— (M — 2), ..., n}, S(n -I- 1) = {n — (M — 3), ..., n -I- 1}. 


If 


max 

S=S+(n-|-l),S_(n-|-l) 


^arg( 2 i) = arg(z,), 

i€S+(n-|-l) 


ieS 


then this maximum cannot exceed , since the sum contains not more then M — 1 summands each 

of which is nonnegative and does not exceed ir/M by our assumption that arg(zi) < tt/M. 

Now let 

max arg(zi) = 2^ -arg(zi). 

i£o 2€S_(n+l) 

Suppose first that n -I- 1 ^ S'_(n -|- 1), i.e., arg(z„+i) > 0. Then S-{n -f 1) C S-{n), implying that 

7r(M- I) 


arg(2i)| < I arg(zi) 

ieiS-{n+l) ieiS-{n) 


< 


M 


by the induction hypothesis. 
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Now consider the case when n + 1 S ^-(n + 1), i.e., arg(z„+i) < 0. We have 

1 


arg(z„+i) = arg an+iw + 


’ ' ' ’ ' ^n—{M—2) 


< 0 . 


Since arg(a„+iw) > 0, we see that 

1 


0 > arg(z„+i) > arg 


If n — (M — 2) G Sj^{n), then 


= - ^ arg(zi) = ^ {-arg(zi)}- ^ arg(z,). 

i=n —(Af—2) i^S—{n) iGS^{n) 


0 > ^ arg(2i) = arg(z„+i) + ^ arg(zi) > - ^ arg(zi) > - 

2 GS_(n+l) iGS-{n) iGS+{n) 

by the induction hypothesis. If n — (M — 2) G S-{n), then 

0 > ^ arg(zi) = arg(z„+i) + ^ arg(zi) - arg(2„_(M-2)) 

iG5_(n+l) iGS—{n) 

7r(M — 1) 


7r(M — I) 
M 


> -arg(z„_(M- 2 )) - arg(z*)>-- 

iGS+ (n) 


M 


□ 


Now suppose the generalized Euclidean algorithm produces all hj positive. Then the sequence 

fn-i+l{z) 

satisfies the conditions of Lemmawhere 0 < arg( 2 :) < tt/M, w = z, ai = hn-i/hn-i+i- Indeed, the 
inequalities (fTlTll hold by Corollary [LTI Each arithmetic polynomial j = 0,... ,M — 1, has 

difference M and degree (M — 1) — j. Hence these polynomials are monomials and therefore all points 


Z{M-l)-j - 


/n —(M—l)+j (^) ^n—{M—l)+j 


/n —(M—l)+j + l (^) ^n—(M —l)+j + l 


for j = 0,1,..., M — 2 lie in the cone If M <i<n then, by the generalized Euclidean algorithm 

and Corollary fTTl 

fn —i ^n—i , fn -i+M{,z) 1 

Zi = --TT = ■;--2 + —- -—=aiZ + 


/n—z+l(^) ^n—i+l fn—i+l{z') 


fn-i + ljz) _ fn-i+2{z) _ ^ ^ fn-i+(M-l){z 
/„_i + 2(z) /„_i+3(z) f„-i + M{z) 


fn-i+l(z) 

fn-i-M(z) 


= QiZ + 


Zi-I ■ Zi-2 ■ ■ ■ Zi_(M-l) 


which implies (I12p . Now using Lemma [521 we obtain the following Corollary. 
Corollary 26. If z G ii^o, 7 r/M and n > M — 1, then 


max y arg , 

S=S+(n),S-{n)\ \fn-i+l{z) 


fn-t{z) 


tt{M -1) 

< -TV- < 

M 


(14) 


The mapping i \—>• n — i maps the set of integers S{n) = {n — {M — 2),..., n} bijectively onto the set 
S{M — 2) — {0 ,..., M — 2}. This transforms (03 into its ’dual’ form 


max arg , 

S=S+(M-2),S_(M-2) I ^ \f^+l{z) 


Mz) 


tt(M — 1) 

< < TT. 

M 


( 15 ) 


Now we can state and prove our main result. 
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Theorem 27. Let f{x) = aoa;" + + • • • + a„ (ao, oi, ..., a„ S R) be a polynomial of degree n 

and M be a positive integer (2 < M < n). Let all the leading coefficients Hq, ..., hn of the polynomials 
/o, ... , fn obtained by applying the generalized Euclidean algorithm with step M, be positive. Then 
f{z) = 0 implies |arg(z)| > 


Proof. Since all coefficients of / are real, its complex roots occur in complex conjugate pairs. Hence 
it is sufficient to prove that f{z) has no roots in the sector 0 < arg( 2 :) < ■^. The equation f{z) = 0 is 
equivalent to 

fo{^) + /i(^) + • • ■ + /m-i(^) = 0, 


or to the equation 


/o 


■(z). 


h 


■{z) + - 


/m-1 /m-1 

We represent each term in the above sum as a product 


fM —2 / \ 

+ 7 - [z) = 

jM-l 


- 1 . 


(16) 


fk fk fk+1 /m—2 

/m-1 fk+1 fk+2 /m-1 


and apply Horner’s algorithm so that the equation (1161) takes the form: 


/m-2 

/m-1 


(^) 


/ /m-3 
\/m-2 




= - 1 . 


(17) 


The proof of the theorem will be completed if we can show that the argument of the left-hand side of CZD 
is smaller than tt in absolute value. As a matter of fact, we will show that it does not exceed tt{M — 1)IM 
in absolute value. To do this we rewrite the left-hand side of m in the form of the following recurrence 
relation 

fk 1 

Uk+i = yk+i{uk + l), ^ k=l,...,M-l, Ui=yi. 

As before, for every k we define the decomposition 


{l,2,...,fc}=r+(A:)UT_(fc), 

where j G T-{k) if and only if arg(yj) < 0 and j G T+{k) if and only if aig{yj) > 0. 
Lemma 28. For every k = 1, ..., M — 1, 

|arg(wfe)|< ^max „ I T" arg(?/,) . 

T=T+(fc),T_(fe) I ^ 


(18) 


Moreover, if arg(ufe) > 0, then 


arg(Mfe) < ^ H 5 

i&T+(k) ^ ’ iGT 


(19) 


if arg(ufc) < 0, then 


-arg(ufc)<- 

iGT-{k) ^ ieT 


( 20 ) 


Proof. We prove the lemma by induction. Since ui = 2/i(0 + 1) it is clear the statement holds for 
k — 1. Suppose that the lemma is proved for k. Then by (TTSl) {yi,y 2 , ■ ■ ■ ,yk,yk+i} is a sequence in 
satisfying 


max 

T=T+(fc+l),T_(fc+l) 


arg(i/,) 

i€T 


Tr{M — 1) 
M 


We have 

arg(ufc+i) = arg(?/fe+i) -|- arg(wfe -|- 1). 
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Suppose first that arg(ufe+i) > 0. Then at least one of arg(j/fe_|_i) and arg(ufc + 1) is non-negative. 
If arg(?/fc+i) > 0 and a,vg{uk + 1) < 0, then 

arg(ufc+i) < arg(yfe+i) < arg(yfe+i) + ^ arg(?/i) = ^ arg(y*). 

If arg(?/fc+i) > 0 and arg(Mfe -|- 1) > 0, then arg(ufc) > 0 and 

arg(ufe+i) = a.rg{yk+i) + arg(ufc -t-1) < a.Yg{yk+i) + a.rg{uk) 

< arg(j/fc+i) -t aigiVi) = Y 

iGT+{k) zGT+(fc+l) 

If aTg{yk+i) < 0 and arg(Mfe -|- 1) > 0, then 

arg(Mfe+i) = arg(j/fc+i) -|- arg(ufc -I- I) < arg(j/fc+i) -|- arg(Mfe) 

< arg(Mfe) < Y arg(zi) = Y arg(yi). 


ieT+(fe) 


i^T+ (fc+ 1 ) 


Suppose now that arg(u/c+i) < 0. Then at least one of arg(j/fc+i) and arg(uk + 1) is negative. If 
arg(2/fe+i) > 0 and arg{uk -I- I) < 0, then 

0 > arg(Mfc+i) > arg(Mfe + 1) > arg('Ufc) > Y arg(j/*) = Y arg(?/,). 

ieT-(k) iGT-{k+l) 

If arg(j/fc+i) < 0 and arg(Mfe -|- 1) > 0, then 

0 > arg(Mfe+i) = arg(yfe+i)-t arg(ufc-h 1) > arg(?/fc+i) > Y arg(yj). 

iGT_(fc+l) 

If arg(?/fc+i) < 0 and arg(Mfe -|- 1) < 0, then arg(ufc + I) > arg(ufe) and 

0>arg(Mfc+i) = arg(j/fc+i)-I-arg(ufc-I-I) > arg(?/fc+i)-I-arg(u/c) 

> arg(j/fe+i)+ Y arg(?/i) = Y 


□ 


jeT_ (fe) 


ieT_(fe+l) 


Returning to the proof of the main theorem, we observe that um-i is the left-hand side of (II3. By 
Lemma [55] and by the inequality (1151) we conclude that um-i cannot equal —1. □ 

Let us illustrate this theorem for n = 5 and M = 3. 

Corollary 29. Let qq, ai, 02 be positive as well as the leading coefficients h^, hi, /15 of the polynomials 
constructed by the generalized Euclidean algorithm with step M = 3 for the polynomial 


f{x) = aox^ -I- aix'^ + a 2 X^ + a^x^ UiX -I- 05 . 


Suppose also that 




oo ai 02 

Then the polynomial f{z) does not vanish in the closed sector 

{^:|arg(^)|<|}. 


( 21 ) 


foix) = aox^ + 032 ;^; 
fi{x) = aix"^ + aix; 
f2{x) = a2X^ + 05 . 
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Proof. First we determine /o, /i, f 2 ' 


Next, 


^ f \ “0 f / \ , ( “004 

fo{x) = —X ■ fi{x) + 03- 

Oi \ Oi 

„ , X Ol f r \ , f “1“5 

fl{x) = —X ■ f2[x) + 04- 

02 V “2 

„ / ^ “102 „ , 

h[X) = - Xfs + 05 

OiOs — 0004 

Then, by (1^ . 

ho = ao > 0, hi = ai > 0, /12 = 02 > 0, > 0, /13 = oq 


fsix) = 1^03 - 
f4{x) = ( 04 - 


(IQCI 4 

Ol 

0105 

02 


fbix) = 05 


03 04 

oo Ol 


> 0 , /14 = Ol 


0-4 O5 
Ol 02 


> 0. 


Notice that m implies that 03, 04 are positive as soon as oq, oi, 02, 05 are positive. 
Using MATLAB or a mere scientific calculator, we can verify that the polynomial 

f{x) = x^ + x'^ + x"^ + l.OOla;^ + a; + 0.999, 

which satisfies the conditions of Corollary [29l has the following roots: 


xi = -1, a:2.3 = -0.49975+ 0.865592i, 0 : 4.5 = 0.49975 ± 0.86617z. 


The roots xi, X 2 , and *3 are located in the left half-plane. The slopes of the vectors 2:4 and X 5 equal 


0.86617 

^0.49975 


±1.73321. 


The slopes of the boundaries of the sector with M = 3 equal 


± tan = ±\/3 = ±1.73205. 

This illustrates numerically the conclusion of Corollary that the polynomial /(z) does not vanish in 
the sector plotted below: 


y 



|arg(z)| < f 


Notice that inequalities m are of course equivalent to inequalities for the following minors of H^: 


ai 

ao 


04 

as 


> 0 , 


02 

Ol 


05 

04 


> 0 . 


TheoremUcan be proved in various ways. The idea of the proof below was suggested by Mikhail Tyaglov. 


Proof of Theorem The proof will proceed by deriving a contradiction from the assumptions 

f{z) = 0, argz € (0, ■^). Since / is a real polynomial, this will also establish by conjugation that the 
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assumptions f{z) = 0, argz € 0) also lead to a contradiction. Finally, / cannot have positive zeros 

since all its coefficients are nonnegative, which will thus show that / has no zeros in the whole sector 
I arg( 0 )| < Disclaimer: this result of course does not prevent / from having zero at the origin! 

Before we embark on the proof, let us make another reduction. Without loss of generality we can 
assume that fo{z) is a polynomial in by multiplying f{z) by a suitable power of z. Of course, this 
may create a (multiple) zero at the origin, but as we already remarked, we are not concerned with those. 
Note that, once fo{z) =:po{z^) is a polynomial in z^, so are z/i(z) =:pi(z^), z^/ 2 (z) =:p 2 {z^), ..., 


,M-1 - 


fM-iiz) ='.pM-i{z^^)- Some of them may be identically zero though. 


Now, let argz G (0, and suppose 


0 = f{z) = fo{z) H-+ fM-i(z) = Poiz‘d) + 


Pi{z 


M\ 


PM-l{z^) 
yM-1 


Z Z" 

Some of the fj, j > 0 may well be identically zero, those can be simply ignored. 

The Hurwitz matrix for each pair {pi,pj), i < j, is a submatrix of the generalized Hurwitz matrix 
Hnif) and is therefore totally nonnegative. Hence each resulting nonzero polynomial pj must have 
only nonpositive real roots by, e.g., [H Theorem 3.44(1)] and Pj/po is an i?-function of positive type, 
i.e., a function mapping the open upper half-plane C+ to itself. A small but crucial technicality: m 
Theorem 3.44(1)] talks about i?-functions of positive type and lists the coefficients of the pair in the 
order opposite to ours here. Since z^ does not lie on the negative real half-line for argz G (0,7r/M), 
neither pj{z^) nor fj{z) can vanish. Furthermore, argz G (0,7r/M), so z^ lies in C-|-, hence 

p(z^) 

arg fj{z) = arg = -j arg z -k </> for some (j) G (0, tt). 

Hence arg/j(z) G {—jn/M, {M — j)7r/M) and fj{z) cannot take on real negative values either. 

As po is certainly nonzero due to its leading term, we can divide everything through by it and get 


Pl(z") P2(.Z^) 


zpoiz^) z^poiz^) 


Pm-i{z^) 

z^-^Po{z^) 


= - 1 . 


( 22 ) 


Among the nonzero terms, some lie in the half-open upper half-plane : ={w : arg re G [0, tt)} and 
some in the open lower half-plane C_ :={w : argw G (—7r,0)}. No terms lie on the negative real axis, 
as we already observed. 

Moreover, the sum of all the terms in is a vector in as well that lies between the positive 
real half-line and the term with the largest argument in the upper half-plane, say, z^~'^pi{z^)/ pq{z^). 
Likewise, the sum of all the terms in C_ is a vector in C_ that lies between the positive real half-line 
and the terms with the smallest (negative) argument, say, z^~^pj{z^^)/po{z^)] j ^ i. This implies that 


arg 


Mzn 


arg 


\zipQ{z^) 


> TT. 


(23) 


and 


Indeed, if not, then all vectors in (ESI lie in the half-plane to the right of the vector Z^ * Pq\z^ 
hence their sum cannot equal —1. 

Again, pi/po, Pjjpo and Pj/pi (if i < j) are i?-functions of positive type by [131 Theorem 3.44(1)], so 

Pi{z^^) 


arg ( 
arg( 


z^Po{z^) 
Pj{^ 


,M\ 


and arg ^ 


^ 


) = —i argz -k (f) 

) = - j arg z -k V' 


for some (^g(0, tt). 
for some ijj € 


z^Poiz^)' 

= (p — Ip, arg = Ip — (p since (p — ip lies within the correct interval (— tt, tt). 


If i < j, then (p — ip G {—tt, 0), hence (j — *) argz + (p — ip e (— tt, {j — i)jp) ■ 

If j < i, then (p-Ip G (O.tt), hence {j - i) argz + (p - ip £ {-{i - - {i- j)w) ■ 

In both cases {j — i) argz + (p — ip cannot exceed tt. But the inequality 

(j — i) arg z + (p — Ip > TT. 

Contradiction! And Theorem]?] is proven. 


means 


□ 
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6 Factorization of the generalized Hurwitz matrix 

The ordinary Hurwitz matrix turns out to admit a particularly simple factorization |12| . with factors 
determined by the leading coefficients arising in the (ordinary) Euclidean algorithm applied to the even 
and odd parts of a polynomial. Our natural question now is whether an analogous factorization holds 
for generalized Hurwitz matrices. This turns out to be true! 

Before we embark on a proof, we will switch from the matrix Hm to its counterpart 




/ ao gm a2M 
0 om -1 a 2 M-i 
0 aM -2 a. 2 M -2 

0 oo gm 
0 0 om-i 


V : 




■7 


This is in keeping with the notation of m-. it also makes for more elegant formulas. The reader intent 
on factoring the original matrix Hm may do so by following the ideas below. 

Theorem 30. Given a polynomial f{x) = gqx'^ + gix^~^ H-|-a„ and a positive integer M >2, let all 

the leading coefficients ho, ..., /i„ of the polynomials /o, ..., fn obtained by applying the generalized 
Euclidean algorithm with step M be nonzero. Then the matrix HM{f) factors as follows: 

HM{f) = J{ci)---J{cn)HM{an), ( 24 ) 

where 

Jic^) = {07“=! > 

{ Ci = ^ ^ if i = j = Mk + 1, fc > 0; 

1 if j = i + 1 , k> 0; 

0 otherwise. 

Proof. We will prove formula (IMl) by induction on n : =deg(/). For deg(/) = 1, we have / = ao(a:)+ai, 
ho = tto, hi = oi, and we obtain 


/ao 

0 

0 

...X 


/an 

ai 

1 

0 •• 

• 0 

0 

...X 


/oi 

0 

0 

...X 

0 

0 

0 



0 

0 

1 • • 

• 0 

0 



0 

0 

0 


0 

Oi 

0 


= 

0 

0 

0 • • 

• 1 

0 



0 

0 

0 


0 

Go 

0 



0 

0 

0 •• 

Qq 

ai 

1 



0 

Oi 

0 


V : 



'■■J 


V : 





■•J 


V : 



•7 


HmU) = J{ci)Hm{gi), 

where ci = — = 7 ^. Hence formula (IM)) holds for n = 1. 
ai ho 

Let us assume that the induction hypothesis holds for n — 1. Represent the polynomial f{x) in the 

usual form f{x) = fo{x) + fi{x) H-h fM-i{x), where each fi picks up the terms with coefficients that 

are i mod M. Run the generalized Euclidean algorithm with step M to generate a sequence /m, • ■ ■, fn- 


21 


















The last polynomial /„ has degree zero. Now define polynomials Fq, ..., as follows: 


Fo{x) = fo{x) + fi{x)-\ -+ fM-i{x){= fix))] 

Fi{x) = fi{x) + f2{x)-\ -h fnix)] 


Fn-{M-l){x) — fn-{M-l){x) + fn-{M-2){x) + • ' ' + fnix)', 
Fn-{M- 2 ){x) = fn-{M-2){x) + fn-{M-Z){x) + • ' ' + fnix)] 


Fn{x) = fn{x). 


Since all the leading coefficients /iq, ..., of the polynomials /o, ..., fn are nonzero, we have deg Fi = 
n — i, i = 0, ..., n. Since /mix) = fo{x) — -r^xfi{x), note that 

hi 


Fi{x) = fi{x) + f2{x) + • • • + (/o(x) - —x/i(x)). 

hi 

We now verify directly that FIm{Fii) = J{ci)FIm{Fi), where ci = — = Since deg(Fi) = n — 1, 

(i\ h\ 

Fi satisfies the induction hypothesis, so the matrix Hm{Fi) factors further as claimed. This proves 
formula (IMl) and the theorem. □ 

Remark. Applying Corollary [121 we obtain that, for j = 1,..., n, 

Hjiiik + 1, r) 


hi = 


HM{k + l,r- 1)’ 


where r = 


M- 1 


Thus, if 

i-l 
M -1 
i — 1 
M -1 


i 


i — 1 

^ hi-i HM{k + 2,r)HMik + l,r-l) 

M - 1 


M -1 

K HM{k+ 2,r-l)HM{k + l,r) 


k = r{M — 1) — i. 


. The case 




M - 1 

is possible if and only if z = l(M — 1) + 1, where I is some nonnegative integer. Then 
i 

— 1, and 


M - 1 


_ HM{k — M + 3,r — 1) ^ F[M{k — M + 3,r — l)HM{k + l,r — 1) 

~ HMik-M + 3,r-2)’ ~ ~ HMik - M + 3,r - 2)HM{k + l,r) 

Just in case of the ordinary Hurwitz matrix, Theorem 1301 yields yet another proof of Theorem 1161 
We restate this theorem slightly to include the matrix HmH) as well. 

Theorem 31. Given a polynomial /(x) = oqx" + aix"“^ + • • • + a„ and a positive integer M > 2, 
let all the leading coefficients ho, ..., hn oi the polynomials /o, ..., fn obtained by applying the 
generalized Euclidean algorithm with step M be positive. Then both matrices HM{f) and HM{f) are 
totally nonnegative. 

Proof. If all coefficients hi are positive, then formula (1241) holds, and all factors in that formula are 
totally nonnegative by inspection. A product of totally nonnegative matrices is totally nonnegative by 
the Cauchy-Binet formula, so the product FIm is totally nonnegative. Since FIm is a submatrix of FIm 
(and, beautifully, vice versa), we conclude that FIm is totally nonnegative as well. □ 


7 Totally nonnegative submatrices of generalized Hurwitz ma¬ 
trices and zero localization 

We now ask another ’natural’ question: is the generalized Hurwitz matrix Hm necessarily totally non¬ 
negative whenever the corresponding (ordinary) Hurwitz matrices h'^m'^ are totally nonnegative for all 
pairs (z, j), 0<z<j<M — 1? The answer to this question happens to be negative. 
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We first recall that the total nonnegativity of the ordinary Hurwitz matrix H 2 {f ) follows from the 
positivity of its leading principal minors up to order n = deg(/) by the classical Hurwitz theorem. We 
will see now that even this (slightly) stronger condition, namely, the positivity of the relevant principal 
minors of all matrices H^\ 0<i<j<M — 1, does not imply that HM{f) is totally nonnegative if 
M ^2. Here is an example. 


Example 32. Let n = 6, M = 3 and f{x) = + 3x^ + 9a;^ + + 2x^ + x + |. Then 

/o(x) = X® + ^x^ + i, /i(x) = 3x® + 2x^, / 2 (x) = Qx"* + X. 

The leading principal minors of the matrix 


= 


/3 2 0 0 

1 2 i 0 

± 29 ^ 

0 3 2 0 

0 1 


V: 


3 1 

2 9 




7 


up to order 4 are 3, 5/2, 4, and 4/9. The leading principal minors of the matrix 


/9 1 0 0 


^( 02 ) _ 


1 I 

0 9 
0 1 


V: 


5 0 
1 0 
3 1 

2 9 


7 


are 9, 25/2, 7/2, and 7/18. Finally, the leading principal minors of 


= 


/9 1 0 • • 
3 2 0 •• 
0 9 1 •• 

Vi i 


up to order 3 are 9, 15, and 15. 
However, 




3 2 0 

1 ^ i 

2 9 

0 9 1 


-2<«- 


Nevertheless, is all matrices have positive leading principal minors for all pairs 0 < i < 

j < M — 1, then all such matrices are totally nonnegative, and this was the only condition (rather than 
the total nonnegativity of the whole matrix) that we used while proving Theorem 0) So, this theorem 
can also be stated as follows: 


Theorem 33. Let /(x) = oox” + oix" ^ + • • • + a„ be a polynomial of degree n and M be a positive 
integer satisfying 2 < M < + 1. If all infinite matrices for i, j = 0, ..., M — 1, i < j, have 

positive leading principal minors, then f{z) cannot have zeros in the cone | arg(z)| < 

We note that Theorem 01 or Theorem [321 also generalize this result [SI Theorem 4.1] for M = n: 

Theorem 34 (Cowling, Thron). Let /(x) = oox" + aix"“^ H-h a„ be a polynomial of degree n with 

all coefficients oq, ..., a„ positive. Then f{z) = 0 implies | arg( 2 ;)| > 
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8 Sufficient conditions for the total positivity of generalized 
Hurwitz matrices 

One last ’natural’ question to address is whether we can prove some kind of a converse to our Main 
Theorem W\ In other words, whether a real polynomial whose zeros lie outside the cone | arg(z)| < 
must have a totally nonnegative generalized Hurwitz matrix Hm- Unfortunately, this happen to be false 
with a vengeance: even the stronger condition that the polynomial / be stable does not imply that the 
generalized Hurwitz matrix Hm must be totally nonnegative for M ^ 2. Here is an example. 

Example 35. Consider the polynomial f{x) = + 5a;^ + 2x^ + 4a: + ^. In this case we have 


H2U) = 


/I 2 
I 5 
0 I 
0 I 
0 0 

V: : 


0 0 

0 0 

5 0 

4 0 

2 - 
^ 2 


7 


The relevant leading principal minors of i? 2 (/) (i-e., the special minors (jb]) for M = 2) are equal to 3, 
5/2, 17/4, and 17/8. Thus, the polynomial f{x) is stable by the Routh-Hurwitz Theorem[TJ 
However, 




A 

2\ 

1 

4 

V2 

V 

1 

2 


= -2 < 0 . 


Thus is not totally nonnegative, and some of the coefficients {/iq, ..., h„} in the generalized 

Euclidean algorithm with M = 3 are negative. 

But in case M is even and / is stable, the total positivity of HM{f) does hold! 

Theorem 36. Let a polynomial/(a:) = aga:"+aia;"’“^ H- |-a„ (oq, oi,..., a„ € R; oq > 0) be stable. 

Tl 

Then, for any M = 2k (k = 1,, — ), its generalized Hurwitz matrix HM^f) is totally nonnegative. 
Proof. Let M = 2k. Consider all special minors 

HmUh) = Hm 

One checks directly that 


j / + 1 
1 2 


j + r - 1 
r 


H. 


M 


j + r-I 
r 


I 




H 2 1 

[k- 

j - 

9 

1 

2 k- 


( 

2 



H 2 1 

[k- 

j - 
2 

2 

2 k- 


j - 1 


J-2 


rk — 


j - 1 


rk — 


2 / 
r + 1 ' 

J-2 


if j is odd; 
if j is even. 


Here A: = y, j = I,..., M - I if r = I, 2,..., 


— 1, J = M — p,..., M — 1 if r = 


M - 1 


(25) 


and 


M -I 

p is the remainder after the division of n by M — I. 

Let us prove that all the minors of 772 (/) defined by Formula (l25ll are positive. Since f{x) is stable 
we have that 772(/) satisfies the conditions of the Routh-Hurwitz Theorem[T]and all its leading principal 
minors are positive. Applying Theorem 3.1 of [5D] or, alternatively. Theorem [T5] fcl. we get 


H. 


*1 


>0 if and only if 0 < 2ji — ii < n, I = 1,... ,r. 


(26) 


Now note that mini<j<r(2(^A: — 


j - 1 


) — 1) = mini<;<j.(/(2fc — 1) — (j — 1)) = M — j if j is odd. 


Likewise, mini<;<r(2(Zfc — 


J-2 


) — (/ -I- 1)) = mini</<r(/(2fc — 1) — (j — 1)) = M — j if j is even. Next, 
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) — 1) = maxi<;<r(^2fc — 1) — (j — 1)) = r{M — 1) — (j — 1) if j is odd. Finally, 
) — (1 + 1)) = maxi<i<r(^2A; — 1) — (j — 1)) = r(M — 1) — (j — 1) if j is even. 


j — 1 

maxi<;<r(2(lfc- — 

j — 2 

maxi<i<r(2(lfc- — 

Since 7=0, ..., M — 1, and M — 1 is odd, we have min (M — j) = 1 and min (M — j) = 2. 
Since 

Tl 

max (r(M- 1) - (j - 1)) = 


M- 1 


(M - 1) - (j - 1), 


we see that 


max( 

j 


n 

M - 1 


(M-l)-(j-l)) 


n 

M - 1 


{M-p-l) 




n 

M -1 


1) +p = n. 


Thus all the minors of the form (1251) are mentioned in formulas (1261) among the principal minors of some 
Hurwitz submatrices above and are therefore guaranteed to be positive. Therefore all special minors of 
Hm^I) are positive, and hence the matrix i?M(/) is totally nonnegative. □ 
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